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INTRODUCTION 


In  seismic  methods  of  yield  estimation  the  amplitudes  of  seismic 
waves  originating  from  an  explosion  must  be  used  to  determine  yields. 

The  two  basic  types  of  seismic  waves  most  frequently  used  are  the 
teleseismic  P  waves  and  the  long  period  Rayleigh  waves  of  intermediate 
(  20  sec)  period.  More  recently  a  considerable  amount  of  work  has  been 
done  to  utilize  the  amplitude  measurements  on  Lg .  For  increasing  the 
accuracy  of  yields  some  estimate  of  the  attenuation  of  the  wave  types  in 
question  is  necessary.  Attenuation  of  20  sec  Rayleigh  waves,  according 
to  the  available  evidence,  does  not  vary  much  from  region  to  region 
since  most  of  the  energy  is  contained  in  the  high  Q  lithosphere, 
although  discontinuities  in  the  lithosphere  may  cause  problems  in  some 
areas.  The  attenuation  o  f  Lg  is  mostly  controlled  by  random 
heterogeneities  in  the  crust  (Aki  1969)  and  will  not  be  discussed  here. 
This  presentation  will  be  devoted  to  the  attenuation  effects  on 
teleseismic,  short  period  body  waves. 

Before  discussing  the  results  of  the  research  concerning  yield 
estimation  from  teleseismic  P  waves  I  want  to  touch  on  some  special 
problems  with  estimating  Q  in  the  short  period  band  that  delayed 
progress  in  this  field  for  a  long  time.  A  more  detailed  discussion  on 
this  topic  by  Cormier  (1982)  can  be  found  in  a  recent  issue  of  the  BSSA. 

The  main  methods  for  measuring  Q  in  the  short  period  band  can  be 
classified  as  follows: 

a)  Relative  spectral  ratio  method,  which  cancels  the  source  effects  for 
axisymmetric  sources  at  various  stations, 

b)  Absolute  modelling  of  the  body  wave  spectra  assuming  some 
"plausible"  source  model." 

c)  Time  domain  modelling  of  the  body  wave  waveforms  using  a  suitable 
source  time  function. 

The  parameters  of  the  path,  the  Q  especially,  must  be  adjusted  to  fit 
the  observations  in  either  the  time  or  frequency  domains  with  a 
frequency  dependent  or  independent  Q  with  some  source  functions.  In  the 


absence  of  noise  all  methods  should  give  the  same  result.  In  practice, 
however,  none  oL  Lliese  conditions  are  met  and  unless  some  provisions  are 
made  to  reduce  the  sensitivity  of  the  method  used  for  the  estimation  of 
Q  to  the  numerous  unknown  parameters  no  meaningful  results  can  be 
obtained.  In  this  presentation  I  shall  emphasize  the  spectral  methods 
at  the  expense  of  time  domain  modelling  of  narrow-band  waveforms  and 
furnish  facts  Lhat  indicate  that  such  methods  are  more  robust,  and  less 
subject  to  bias  due  to  the  various  still  unresolved  basic  problems 
concerning  seismic  source  and  site  related  distortions  of  the  spectra 
and  the  lack  of  the  precise  knowledge  about  the  possible  frequency 
dependence  of  Q. 

In  Slide  1  some  theoretical  shapes  of  P  wave  spectra  are  shown 
computed  for  a  200  kt  explosion  source  using  the  VSB  source  model.  The 
sensitivity  of  the  shapes  to  t*  is  quite  evident  and  small  differences 
in  t*  are  quite  discernible  in  the  t*  range  of  0-.3.  At  higher  t* 
values  tho  accuracy  of  t*  measurement  must  decrease  because  of  the 
limited  frequency  range  available  for  the  fit  and  the  sensitivity  of  the 
details  in  the  spectra  to  other  perturbing  factors. 

The  key  to  the  estimation  of  the  attenuation  is  the  extreme 
sensitivity  of  the  spectral  shapes  to  Q,  at  the  expense  of  other 
factors.  Assuming  losses  in  shear  deformation  in  the  mantle,  which 
appears  to  agree  with  most  experimental  facts,  and  using  some  basic 
physical  limits  on  plausible  source  functions  it  is  possible,  using 
spectral  analyses,  to  put  some  limits  on  the  possible  range  of  Q  and  its 
frequency  dependence. 

To  illustrate  the  insensitivity  of  short  period  P  wave  spectra  to 
near-site  focusing  effects  one  can  compare  the  variability  of  the  P  wave 
amplitudes,  waveforms  and  shades  of  spectra  at  seismic  arrays.  While 
seismic  waveforms  and  amplitudes  exhibit  an  extreme  variation  from 
sensor  to  sensor  (across  NORSAR  variations  of  5:1  in  amplitudes  are 
common)  estimates  of  t*  from  common  sources  at  the  various  sensors  have 
a  standard  deviation  of  the  order  of  .06  sec  if  one  considers  the  .5-4 
Hz  band  at  the  same  array.  (Slides  2  and  3)  It  would  not  be  wise, 
therefore,  to  rely  on  absolute  P  wave  amplitude  variations  to  estimate 


relative  Q.  The  amplitude  variations  in  short  period  P  waves  have  been 
shown  to  be  the  results  of  focusing  on  lateral  heterogeneities  in  the 
crust  and  upper  mantle  under  the  array  and  is  thus  not  related  to  Q. 
(Slide  4)  This  reasoning  applies  to  all  measurements  of  Q  at  individual 
recording  sites  using  wave  amplitudes.  Only  averages  of  amplitudes  from 
many  azimuths  may  be  meaningful. 

Another  natural  consequence  of  small  scale  lateral  heterogeneities 
in  the  crust  is  that  the  seismic  rays  are  diverted  near  the  source  to 
constitute  large  scale  focus ing-de focus ing  patterns  at  teleseismic 
distances  (Hadley  1979,  Butler  and  Ruff  1980).  An  example  of  such 
patterns  after  Hadley  (1979)  is  shown  in  Slide  5. 

Similarly  the  possible  range  of  Q  and  its  frequency  dependence  may 
be  delimited  by  the  range  of  observable  signal  frequencies  and  by  the 
shapes  of  spectral  ratios  if  frequencies  up  to  10  Hz  are  used.  It  can 

be  shown,  for  instance  that  the  spectra  and  waveforms  of  short  period  S 

waves  put  some  quite  tight  upper  limits  on  t*  along  some  paths,  even 
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assuming  a  delta  function  source.  The  shapes  of  wide  band  spectra  in 
the  short  period  band  do  not  seem  to  support  any  rapid  variation  of  Q 
with  frequency,  but  frequency  dependence  of  Q  is  necessary  to  reconcile 
long  period  Q  estimates,  if  they  are  to  be  believed,  with  those  obtained 
in  the  short  period  band.  (Slide  6a,  b) 

In  general  studies  of  Q  in  narrower  bands  do  not  support  rapid 
variations  of  Q  with  frequency  (Der  et  al  1980,  Lay  and  Helmberger  1981. 
Shore  1983).  The  arguments  for  frequency  dependence  of  Q  are  more 
convincing,  however,  if  one  desires  to  reconcile  long  and  short  period  Q 
estimates,  but  even  in  this  case  Q  is  very  poorly  constrained  in  the 
long  period  band  for  specific  ^regions .  Estimates  from  free  oscillation 
may  also  be  seriously  biased  relative  to  the  real  Q  structures  in  the 
most  typical  regions  of  the  world  because  of  severe  lateral  variations 
in  Q  and  velocities.  Therefore  the  question  of  frequency  dependence  of 
Q  can  not  be  decided  until  more  research  is  done,  but  it  seems  certain 
that  any  changes  of  Q  with  frequency  is  quite  gradual  along  most  paths. 


In  view  of  the  absence  of  indications  of  strong  frequency 
dependence  1  find  Lhe  attempts  of  Lundquist  and  Samowitz  (1982)  to  fit 
detailed  pairs  of  absorption  band  models  to  interstation  P  wave  spectral 
ratios  quite  '  renin  t  n  r  e  ,  and  only  careful  ,  broad  band  studies  of  Q 
extending  Lo  low  frequencies  will  be  able  to  determine  the  details  of 
frequency  dependence,  instead  of  the  narrow  band  data  used  by  Lundquist 
and  Samowitz  (typically  .5-4  z.). 

If  most  of  the  anelastic  losses  occur  in  the  upper  mantle  t*  should 
not  be  strongly  dependent  on  distance  along  paths  in  laterally 
homogeneous  structures.  This  appears  to  be  the  case  as  long  as  we 
attempt  to  keep  the  structure  separated.  Shore's  (1983)  results  and  our 
own  show  the  weak  dependence  of  t*  on  distance.  The  somewhat  higher 
values  for  these  t*  for  shields  are  due  to  the  fact  that  they  were  tied 
to  PP  measurements  which  because  of  the  imperfect  surface  reflection  in 
the  short  period  band  are  likely  to  give  upper  limits. 

For  nuclear  explosions  the  presence  of  surface  refections  and 
possible  spall  arrivals  are  a  little  more  than  a  slight  perturbation  in 
the  estimated  Q  value  as  long  as  one  assumes  a  rather  weak  scaling 
effect  on  the  spectra  and  a  wide  band  of  frequencies  are  considered, 
and  no  matter  what  the  pP  parameters  are,  the  low  and  high  t*  paths 
yield  spectra  that  are  quite  distinguishable. 

In  view  of  the  extreme  site  variations  in  the  short  period  band  it 
is  obvious  ilia  i  absolute  amplitude  modelling  of  waveforms  is  much  less 
accurate  for  modelling  Q.  Likewise,  the  determination  of  other  source 
parameters  using  absolute  waveforms  has  been  found  quite  nonunique 
Cormier  (1982). 

Although  many  papers  in  t^he  literature  that  use  t*  *  1.0  routinely 
in  the  short  period  band,  regardless  of  the  tectonic  setting  of  the 
path,  such  results  are  usually  quite  poorly  constrained  (Der  and 
McElfresh  1980).  Moreover  a  lower  t*  would  usually  fit  the  data  as 
well.  In  most  cases  t*  trades  off  with  source  parameters,  In  one  such 
study  (Burdick  and  Helmberger  1977)  the  derived  overshoot  ratios  were 
appropriate  to  a  source  medium  of  the  consistency  of  rubber. 


Exclusively  high  values  ( t*  *  1.15)  by  Helraberger  and  Handley  would  lead 

to  unreasonably  high  bias  estimates. 

As  a  general  comment  I  wish  only  to  state  that  near  and  above  1  Hz 

the  value  of  t*  =  1  appears  to  be  quite  inappropriate  for  most  paths, 

not  only  shields.  Results  for  shields  give  quite  low  absolute  t*  from 

PP  (Shore  1983)  or  t*  from  short  period  S  phases  (Der  et  al  1982). 

Similar  results  were  obtained  by  Stewart  (1983)  from  rise  times  of  P 

waves  ( t*  <  . 5  sec )  . 

P 

In  the  following  discussion  of  the  effect  of  mantle  Q  on  yield 
estimation  we  shall  make  the  following  implicit  assumptions. 

a)  The  frequency  dependence  of  Q  is  weak  along  the  various  paths. 
Therefore  the  Q  related  amplitude  variations  will  be  correlated  with 

the  variations  of  the  shapes  of  P  wave  spectra. 

b)  We  assume  that  using  large  number  of  relative  spectral  ratios  among 
recordings  at  various  stations  the  source  effects  will  cancel. 

c)  That  explosions  in  hard  rock  will  have  source  spectra  adequately 
described  by  the  von  Seggern-Blandford  model. 

These  assumptions  agree  with  most  observed  facts  we  know  of  at  the 
present  time.  Exceptions  appear  to  occur  for  explosions  in  salt,  clay, 
and,  for  high  frequencies,  at  Yucca  Flats,  Within  the  above  constraints 
we  shall  attempt  to  give  our  best  estimates  of  the  differences  in  the 
mantle  attenuation  under  several  major  test  sites  and  describe  any  data, 
less  complete,  for  others.  The  estimates  will  be  based  on  the  best  fits 
to  available  estimates  of  absolute  and  relative  apparent  t*. 

After  these  preliminaries  I  shall  review  now  the  available  evidence 
indicating  major  differences  of  mantle  Q  under  some  of  the  most 
important  test  sites. 

The  problem  of  primary  concern  in  yield  estimation  is  the  so  called 
"magnitude  bias",  i.e.,  variations  of  the  m^-yield  relationships  for 
various  test  sites  caused  by  differing  amount  of  amplitude  reduction  of 
P  waves  under  various  test  sites  due  to  mantle  attenuation. 


Another  problem  of  great  interest  in  nuclear  monitoring  is  that  of 
U  i  sc  r  i  in  i.  na  i.  i  u  u  .  Whiio  tor  explosions  oi  largo  yield  the  M^-iu^ 
discriminant  is  very  effective,  at  lower  yields  the  seismic  surface 
waves  may  mu  l>o  detectable.  To  d  i  sc  r  i  in  i  na  t  e  at  such  yields  the  use  of 

the  relative  spectral  contents  of  explosions  and  earthquakes  was 
proposed.  Explosions,  being  sources  of  small  dimensions,  can  be 
expected  to  generate  low  frequency  waves  less  effectively  than 
earthquakes.  Unfortunately,  the  spectral  shapes  are  also  affected  by  Q 

and  the  a  priori  knowledge  of  Q  along  the  path  is  necessary  for 
effective  discrimination. 

The  mantle  Q  along  the  path  is  also  a  major  factor  in  the 
information  content  of  teleseismic  signals.  The  accuracy  in  determining 
the  details  of  the  source  functions  or  the  amplitudes  and  delay  times  of 
surface  reflections  may  be  significantly  reduced  if  the  available  signal 
bandwidth  is  diminished  by  low  Q  along  the  path.  This  will,  in  turn, 
affect  the  reliability  of  any  correction  of  the  amplitude  for  the 
surface  reflection  and  thus  the  estimated  yield  itself. 

Seismological  evidence  that  Q  variations  pose  an  important  problem 
in  yield  estimation  emerged  fairly  early  in  the  studies  of  seismic  waves 
from  U.  S.  Nuclear  explosions.  Seismic  waves  from  the  nuclear 
explosions  Gnome  and  Salmon  showed  a  strange  asymmetry  in  the  amplitudes 
and  frequency  contents  of  the  observed  P  waves  (Slide  7)  which  was 
noted  in  the  sixties  but  largely  forgotten  or  ascribed  to  differences  in 
crustal  structures  soon  after.  Many  recent  workers  agree  that  the  first 
arrivals  must  have  travelled  through  the  upper  mantle  beyond  a  certain 
distance  from  these  events  so  that  the  reduction  in  the  high  frequency 
contents  of  P  in  the  western  United  States  but  not  in  eastern  North 
America  reveals  a  difference  in  mantle  attenuation.  Plots  of  apparent 
t*  values  derived  from  the  estimated  source  spectrs  of  these  events  show 
a  constant  slope  characteristic  of  attenuation. 

P  wave  spectra  of  explosions  at  the  Nevada  Test  Site  were  shown  in 
early  studies  to  be  characterised  by  a  general  lack  of  high  frequency 
content  beyond  4-5  Hz  at  teleseismic  distances  overall,  while  most 
explosions  in  the  shield  areas  of  Eurasia  and  recorded  in  the  shield, 


were  shown  to  have  significant  high  frequency  energy  in  their  P  waves 
with  8-10  11  v.  energy  still  above  the  noise  level. 

The  path  related  spectral  anomalies  in  P  waves  were  also  associated 
with  anomalies  in  the  amplitudes  of  the  observed  teleseismic  P  waves 
across  North  America  reported  in  several  major  studies.  The  results  of 
Booth,  Marshall  and  Young  (1974)  demonstrated  these  for  P  arrivals.  An 
apparently  conflicting  conclusion  was  made  by  Butler  and  Ruff  (1980) 
asserting  that  there  are  no  m^  differentials  between  the  m^  values 
observed  in  the  EUS  and  WUS  or  the  US  as  a  whole.  Proper  consideration 
of  the  total  pattern  of  magnitude  anomalies  across  the  U.S.  (Douglas  and 
Marshall  1983)  however  totally  invalidates  this  claim.  It  so  happens 
that  Butler  and  Ruff  considered  the  averages  of  the  west 
coast-intermountain  belt  and  the  NEUS-EUS  regions  respectively  in  which 
the  first  parL  was  always  over  represented  in  station  coverage.  I  think 
nobody  after  considering  Slide  8  which  shows  the  ra^  results  of  the 
merged  LRSM-WWSSN  stations,  could  rightly  claim  that  no  m^  anomalies 
exist  across  the  United  States.  Studies  of  S  wave  arrivals  showed 
similar  regional  patterns  of  amplitude  and  spectral  anomalies  but  much 
accentuated  in  a  manner  suggestive  of  energy  losses  in  shear  deformation 
in  the  mantle.  Attempts  to  model  these  amplitude  and  spectral  anomalies 
with  variations  of  crustal  effects  failed.  The  amplitude  anomalies 
could  not  be  explained  by  the  effects  of  even  extreme  crustal  models 
(Slide  9).  Note  in  this  figure  that  the  average  WUS  -  EUS  m^  difference 
is  0.33  and  that  RKON  is  in  the  middle  of  the  EUS  population  -  not 
the  highest  amplitude  station  by  any  means. 

Another  early  indication  that  the  m^-yield  scales  depend  on  the 
source  region  comes  from  studies  of  the  M^m^  relationships  for 
explosions  (Slide  10)  which  showed  a  too  low  m^  for  NTS  explosions  (or 
alternately  Loo  high  M^).  This  can  easily  be  explained  by  greater 
attenuation  under  the  WUS  of  the  P  waves.  Simulations  confirmed  the 
regression  results  (Slide  11)  and  showed  that  a  formula  relating  ra^  to 
t*  Amb  *  1.35  At*  is  approximately  valid.  The  formula  was  empirically 
derived  by  fitting  regressions  to  WUS  and  EUS  crustally  corrected  A  m^ 
values  by  Der  et  al.,  (1979). 


Since  in  estimating  yields  we  are  more  interested  in  the 

geophysical,  properties  of  a  few  test  sites  rather  than  generalities 

about  regional  variations,  during  the  late  seventies  ARPA  undertook  a 

major  effort  to  measure  the  size  of  the  intersite  "magnitude  bias" 

between  NTS  and  a  hypothetical  shield  station.  Since  the  Soviet  test 

sites  are  not  accessible  to  us  the  shield  site  was  typified  by  the 

station  RKON  (Red  Lake,  Ontario)  on  the  southern  edge  of  the  Canadian 

shield.  The  stations  HNME  and  IFME  in  Maine  were  found  to  be  similar  in 

surface  geology  to  the  Kazakh  test  and  were  therefore  selected  as 

possible  analogs.  These  sites  were  also  occupied  at  various  stages  of 

the  experiment.  At  NTS  we  had  two  stations  on  the  Climax  stock,  one  at 

the  Gold  Meadows  stock,  two  at  Pahute  Mesa  and  four  at  Yucca  Flats. 

Later  more  stations  located  at  various  nuclear  test  sites  in  the  western 

United  States  were  added.  These  occupied  the  sites  of  the  Gasbuggy,  Rio 

Blanco  and  the  Faultless  explosions.  The  site  of  another  explosion  in 

granite,  Shoal,  was  occupied  earlier  by  an  LRSM  station.  (Slide  12) 

After  approximate  corrections  for  the  crustal  structure  were  applied  the 

trace  amplitudes  revealed  lower  amplitudes  at  all  UUS  sites  relative  to 

RKON.  (Slide  13)  The  m^  differences  (log  A/T)  are  reduced  from  these 

single  (log  A)  differences.  The  measurements  of  apparent  t*  are 

especially  revealing  (Slide  14).  There  appears  to  be  a  significant 

difference  in  the  frequency  contents  of  all  teleseism ic  arrivals  at  the 

WUS  sites  relative  to  RKON.  The  two  sites  in  the  NEUS,  HNME  and  IFME, 

appear  to  have  an  intermediate  position  between  the  shield  and  the  WUS. 

As  we  shall  see  the  Kazakh  t*  values  are  low  so  that  these  intermediate 

t*  values  suggest  that  the  NEUS  stations  HNME  and  IFME  are  not  good 

analogs  of  the  Kazakh  test  site.  The  differences  in  the  frequency 

contents  of  P  waves  are  also  apparent  in  the  average  differentials  in 

dominant  period  at  stations  with  similar  instrument  responses. 

* 

Similar  regional  dependence  in  t*  was  also  found  for  P  wave 
arrivals  from  Soviet  nuclear  explosions  across  North  America  (Slide  15) 
by  Shore  (1983).  For  this  data  set  the  RKON-Basin  and  Range  t* 
differential  is  greater  than  in  the  SDCS  experiment. 

This  experiment  has  found  that  the  test  sites  within  the  WUS  fit 
well  into  the  regional  pattern  of  mantle  attenuation  in  the  same  region. 


it 
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The  average  t differential  between  RKON  and  the  station  at  the  Climax 
Stock  ;il  MTS  was  near  .2  sec  which  corresponds  to  about  a  0.27  magnitude 
differential  assuming  a  constant  Q. 

The  patterns  of  regional  attenuation  found  in  the  NTS  experiment 
were  furLher  confirmed  by  and  fit  well  into  the  results  of  broader 
regional  studies  of  mantle  Q  under  the  United  States  (Der  et  al.,  1982, 
Lay  and  Helmberger  1980).  The  studies  of  short  period  S  waves  yielded 
especially  dramatic  examples  that  confirm  the  reality  of  lateral  Q 
variations.  There  can  be,  therefore,  hardly  any  doubt  about  the 
existence  of  lateral  variations  of  Q  under  the  United  States.  Slide  16 
shows  a  good  example  of  variations  in  the  amplitudes  and  wave  periods  of 
short  period  S  waves  from  deep  earthquakes  across  the  United  States  due 
to  lateral  changes  in  mantle  Q.  Burdick  (1982)  has  laso  concluded  that 
most  of  the  United  States  is  underlain  by  one  of  two  types  of  mantle. 

Following  the  NTS  experiment  we  have  undertaken  an  extensive  study 
of  the  spectra  of  Leleseismic  P  waves  from  NTS  and  Kazakh  nuclear 
explosions  as  recorded  at  a  large  number  SRO,  LRSM  and  array  stations. 

We  are  now  dealing  wiih  a  different,  more  difficult,  problem; 
determination  of  t*  from  the  spectra  by  assuming  the  source  spectrum 
instead  of  cancelling  the  source  relative  measurements. 

The  von  Seggern-Blandford  model  was  used  as  a  preliminary  source 
model.  Despite  the  possiblity  of  multiple  arrivals  and  other 
complications,  comparison  of  the  spectra  from  explosions  with  varying 
estimated  yields  seemed  to  obey  the  cube-root  scaling  law  we  have  used 
for  the  von  Seggern-Blandford  model.  Any  uncertainties  of  the  assumed 
source  specLra  will,  in  this  case,  affect  the  t*  estimates.  We  have 
accumulated  more  than  100  spectra  during  this  experiment.  By  comparing 
the  spectra  of  NTS  and  Ka zakh 'exp  1  os i ons  for  events  with  comparable 
yields  the  major  difference  in  the  spectral  contents  in  P  waves  from  the 
two  test  areas  is  quite  obvious.  (Slides  17a,  b  show  examples  from 
NORSAR) .  At  arrays  the  overall  spectral  shapes  remain  consistent  large 
despite  site  amplitude  variational 


In  order  to  interpret  the  total  data  set  and  to  exploit  any 
iutciiwil  cons  i  sLuuc  ius  in  tlio  data  we  liave  assumed  the  simple 
statistical  model;  _  *•*  +  t*  + 


d*. 

i.l 


c? 


eij 


Where  tjj  is  the  measured  average  t*  from  test  site  i  at  station  j,  tij 

and  t*  are  the  test  site  and  station  terms  and  e. .  is  an  error  term.  In 

i  i  J 

this  model  we  ignore  the  distance  dependence,  which  at  teleseismic 

distance  is  small,  and  any  dependence  on  azimuths.  We  therefore  attempt 

to  determine  average  upper  mantle  Q  contributions  to  t*  under  the  test 

sites  and  receivers  averaged  in  the  small  cone  in  the  mantle  sampled  by 

teleseismic  ray  paths. 

Since  the  test  site  and  station  terms  trade  off  against  each  other, 
no  unique  solution  exists  unless  we  impose  some  constraint.  Since  in 
the  SDCS  experiment  we  had  a  measured  t*  difference  between  Pahute  Mesa 
and  RKON  of  .2  sec  we  have  adjusted  the  differences  in  the  regression 
results  to  this  independently  determined  value.  (Actually  this 
differential  was  measured  more  accurately  than  the  accuracy  of  the 
regression  results). 

By  looking  at  spectra  of  Yucca  explosions  at  regional  distances  we 
have  concluded  that  the  von  Seggern~Blandford  granite  RDP  used  in  our 
analyses  is  not  valid  for  these  events.  We  decided  therefore  to  use  the 
Yucca  data  only  for  refining  the  station  terms  (which  are  determined  by 
spectral  ratios).  Determination  of  Yucca  RDP-s  requires  more  work. 
Actually  the  direct  Yucca  t*  estimates  are  not  needed  to  determine  Yucca 
test  site  t*  because  we  have  a  very  large  data  set  of  reciprocal 
measurements  that  show  that  Yucca  Flats  are  not  significantly  different 
from  Pahute  Mesa  with  regard  to  mantle  Q.  The  regression  on  the  other 
hand  shows  a  rather  large  At*#between  Pahute  Mesa  and  the  Kazakh  test 
sites.  (Slide  18) 

The  station  terms  on  the  left  do  not  vary  much,  although  WHYK  and 
ZOBO  are  high.  NORSAR  and  NPNT  are  low.  The  test  site  terms  show  a 
great  difference  between  NTS  and  Kazakh  (Degelen  and  Shagan) .  Based  on 
these  results  one  can  say  with  certainty  that  NTS  and  Kazakh  must  have 
quite  sizeable  differences  in  the  attenuative  properties  of  the  upper 


mantle,  but  there  may  be  a  sizeable  uncertainty  in  the  NTS  measurements 
due  to  the  narrow  band  P  wave  spectra  from  NTS.  The  station  terms  are 
less  affected  by  this.  Therefore  we  cannot  base  our  estimate  of  the 
bias  between  the  two  test  sites  on  the  results  of  these  direct,  not 
relative,  t*  estimates  alone. 

The  regression  gives  a  higher  At*  bias  for  Pahute  Mesa  relative  to 
Kazakh  than  the  0B2NV-RK0N  difference.  Can  it  be  that  Kazakh  and  NORSAR 
have  higher  Q  in  the  underlying  mantle  than  RKON?  This  appears  to  be 
the  case  if  one  compares  the  spectra  of  P  waves  from  common  events  at 
RKON  and  NORSAR.  The  corresponding  differential  in  t*  is  about  .1  sec. 
Therefore  it  is  possible  to  substantiate  a  NTS-Kazakh  bias  of  the  order 
of  .3  in  t*.  In  the  course  of  making  t*  estimates  from  Pahute  Mesa  NTS 
explosions  we  did  not  note  any  major  deviation  in  the  source  spectral 
shape  from  the  von  Seggern-Blandford  model,  and  thus  the  .3  bias  seems 
to  be  reasonable  considering  all  the  evidence.  Another  relevant 
question;  can  it  be  that  the  Kazakh  explosions  put  out  anomalously  high 
amounts  of  high  frequency  energy?  This  however,  can  be  ruled  out 
outright  by  considering  the  extremely  high  frequencies  in  the  P  wave 
spectra  up  to  10  llz  that  would  be  impossible  to  observe  with 
substantially  higher  t*  even  if  one  assumed  a  flat  source  spectrum.  The 
t*  estimates  from  Soviet  PNE-s  are  essentially  similar  to  those  obtained 
from  Kazakh  explosions. 

An  equivalent  of  the  reciprocal  experiment  performed  at  NTS, 
although  not  as  exhaustive,  was  made  possible  by  published  data  from  a 
Soviet  CISS  station  near  Semipalatinsk .  Comparison  of  the  spectra  of 
Alaskan  earthquakes  at  NTS  and  at  this  station  resulted  in  a  bias 
estimate  in  t*  of  .24  sec  by  Murphy  and  Tzeng  (1982). 

Blandford  (1981)  comparing  carefully  selected  Soviet  explosions  to 
Piledriver  at  the  common  station  NPNT  obtained  a  At*  of  0.25  as  a  best 
estimate . 

Burdick  (1982)  comparing  the  WUS  and  the  Eurasian  shield  via  S-SS 
long  period  delays  concludes  that  the  t*  difference  between  these  broad 
areas  should  be  0.2. 


Considering  the  total  evidence  it  can  be  said  that  the  NTS-Kazakh 
t*  bias  duo  Lu  excessive  attenuation  under  the  basin  &  Range  must  be  in 
the  range  of  .2-. 3  sec.  It  is  probable  that  it  does  not  exceed  .3  but 
is  almost  certain  to  be  more  than  .2.  Assuming  the  formula  *  1.35 

At*  the  magnitude  difference  corresponding  to  this  t*  difference  would  be 
0.27  to  0.40  in  m^ .  However,  in  the  SDCS  experiment  we  noted  that  the 
m^  differential  between  NTS  and  RKON  was  0.1  less  than  that  implied  by 
the  At*  values.  This  was  attributed  to  focusing  under  NTS.  This 
correction  must  also  be  applied  in  the  present  study  so  that  our  final 
bottom  line  is  a  Am^  in  the  range  0.17  to  0.30  m^ . 


STHDT ES  OF  P  WAVES  FROM  OTHER  TEST  SITES 


Algeria 

A  site  of  considerable  interest  in  yield  determination  is  the 
Hoggar  massiv  in  Algeria  where  some  of  the  French  nuclear  tests  took 
place.  This  site  is  presumably  locaLed  over  a  so  called  "hot  spot"  and 
i 8  characterized  by  incipient  rifting,  recent  volcanism  and  thinning  of 
the  crust.  Comparison  of  P  wave  spectra  from  the  nuclear  explosion 
Rubis  at  Hoggar  and  a  Kazakh  explosion  at  the  AWRE  array  EKA  indeed 
shows  a  decrease  in  the  high  frequency  content  of  P  waves  from  Hoggar 
compared  to  those  from  Kazakh.  This  supports  the  idea  that  Hoggar  is 
indeed  located  over  a  "hot  spot".  Observations  of  P  waves  spectra  of 
nuclear  explosions  at  Hoggar  also  show  the  same  decrease  of  high 
frequency  content  in  P  waves  relative  to  shield  type  of  paths. 

Comparing  the  apparent  t*  of  the  nuclear  tests  at  Hoggar  to  those 
at  NTS,  the  Hoggar  site  appears  to  be  intermediate  between  the  shield 
and  NTS  with  respect  to  attenuation  in  the  underlying  mantle. 

Novaya  Zemlya 

Novaya  Zemlya  nuclear  explosions  observed  at  high  Q  stations  yield 
low  t*  estimates  indicating  that  the  Q  in  the  mantle  under  this  source 
region  is  high,  and  comparable  to  that  under  Kazakh. 

Lake  Baykal  Area 

Seismic  P  waves  recorded  at  LASA  and  NORSAR  for  two  nuclear 
explosions  in  the  Baykal  area  did  not  show  any  indication  of 
encountering  a  low  Q  upper  mantle  along  these  paths.  (Sobel  et  al., 
1977)  .  In  a  more  detailed  st^dy  covering  a  wider  range  of  azimuths 
Savino  et  al  .  ,  (  1975  )  have  found  that  the  spectral  discriminant  failed 
along  some  azimuths  (Slide  19).  This  can  be  interpreted  as  the  effect 
of  low  mantle  Q  under  the  rift  zone,  also  indicated  by  other  geophysical 
studies . 


I  ml  i.  a  n  N  u  a  !  a  a  i 


M x plosio n 

The  rise  times  of  P  waves  and  spectral  analyses  of  P  from  the 

Indian  u  ml  a  r rou  ml  explosion  at  NORSAR  are  consistent  with  a  high  Q 
path  . 


Soviet  PNE-s 


As  we  have  already  mentioned  Soviet  PNE-s  must  be  located  in  the 
areas  of  high  mantle  Q  considering  the  rise  times  of  P  and  their 
spectra.  This  includes  known  shots  in  salt. 

French  Nuclear  Explosions  in  the  Pacific 

These  have  not  been  studied  much  with  respect  to  Q  but  published  P 
wave  spectra  at  the  Hagfors  array  in  Sweden  (Nedgard,  1978)  are 
indicative  of  high  mantle  Q  under  these  test  sites. 

Amchitka  Tests 

Analyses  of  the  U.  S.  nuclear  tests  on  Amchitka  island  usually 
result  in  t*  at  stations  in  the  contiguous  U.  S.  of  the  order  of  .4— ,5. 
Considering  the  fact  that  the  mantle  Q  under  some  U.  S.  stations  is  low, 
this  indicated  low  or  moderate  attenuation  in  the  source  area.  These 
admittedly  scarce  results  agree  with  the  findings  of  Barazangi  et  al . , 
(1975)  that  there  is  no  extensive  back-arc  attenuation  zone  in  the 
mantle  near  the  Aleutian  islands. 
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CONCLUSIONS 


/ 


Spectral  analysis  of  P  waves  crossing  the  mantle  under  various  test 
situs  has  rovu.tlotl  considerable  variations  in  the  altenu.it  ive  properties 
of  the  mantle  which  must  be  considered  when  computing  yields  of  nuclear 
explosions.  The  t*  bias  between  the  NTS  and  Kazakh  test  sites  must  be 
in  the  lange  0.2  to  .3  sec.  There  are  indications  of  anomalous  zones  of 
attenuation  under  the  French  test  sites  in  the  Sahara  and  near  some 
suspected  nuclear  explosions  in  the  Baykal  area.  ^The^jjjant le  attenuation 
under  the  Sahara  test  sites  appear*  to  be  intermediate  between  a  shield 
and  the  western  United  States.  <:^::Fh3>  other  test  sites  for  which  we  have 
relevant  results  seem  to  be  in  areas  of  low  mantle  attenuation. 

The  patterns  of  attenuation  found  by  the-analyses  of  nuclear  explosion 
data  fit  roughly  in  the  worldwide  patterns  found  in  more  general  studies 
of  plate  tectonics.  No£e--ef  they  test  sites  are  located  in  areas  of 
extremely  low  Q  zones  in  the  mantle,  such  as  those  found  under  some 
back-arc  basins. 


Empirical  analyses  and  synthetic  simulations  give  the  formula  J 

no. 

|,3  r  Q>-'  S  t*C*U  Is  p) 


"b  ' 


1.35 


A  t* 
P 


for  deriving  the  m,  bias  from  a  t*  differential. 

'P 

r  (yS-Vfc,  Ic)  j 

“"The- regional  patterns  of  lateral  Q  variations  found  by  spectral 

analyses  are  'a-lso  -well  correlated^wx  th  «-e  mmrher  ef -other  relevant 
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geophysical  parameters,  such  as  P.  velocities,  travel  time  variations  of 

A 


P  and  S  waves,  heat  flow,  and  upper  mantle  conductivity . ^Although  care 
must  be  exercised  in  the  use  of  such  quantities  for  deriving  mantle  Q 
since  these  other  parameters  are  not  directly  related  to  Q,  they  can  be 
used  as  diagnostics  in  the  absence  of  direct  data. 


Indirect  methods  of  estimating  the  mantle  attenuation  using  body 


wave  (?n,  S  and  SS  velocities  have  been  proposed  (Marshall,  Springer  and 


Rodean  1979,  Burdick  et  al.,  1982)  and  may  be  of  some  value  in 
inaccessible  regions.  (Slides  20  and  21).  Strictly  speaking,  however, 
body  wave  velocities  may  only  be  loosely  related  to  Q  and  there  are 
indications  that  body  wave  travel  times  worldwide  may  not  correlate  with 
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nuclear  explosion  with  t*.  At  low  t*  these  shapes  are  very 
diagnostic,  but  the  accuracy  of  estimation  decreases  with 

increasing  t*. 

2  Illustration  of  the  variations  of  P  wave  amplitudes  across 

two  arrays. 

3  Scatter  in  apparent  t*  across  NORSAR.  While  the  P  wave 
amplitude  commonly  vary  by  a  factor  of  five,  the  apparent  t* 
have  a  low  (.06  sec)  standard  deviation. 

4  Azimuth-distance  dependent  P  wave  focusing  patterns  at  LASA. 
After  Chang  and  von  Seggern  (1981). 

5  Azimuthal  variations  in  amplitudes  for  various  NTS  events 
indicating  near  source  focusing  effects  (After  Hadley  1979). 

6a, b  Some  t*-vs- frequency  curves  (a)  and  the  associated  spectral 

ratio  shapes  (b) .  Small  changes  results  in  drastic  and 
detectable  changes  in  the  spectral  ratios  putting  severe 
limits  on  acceptable  t*-vs-frequency  relationships.  In 

genera;  rapid  changes  can  be  ruled  out. 

7  P  waveforms  for  the  Salmon  nuclear  explosion.  Variations  in 
the  frequency  content  at  upper  mantle  distances  reflect 
lateral  variations  in  Q. 

8  Station  terras  for  M.  for  the  merged  LRSM-WWSSN  networks 
across  the  United  States  (Douglas  and  Marshall  1983). 

9  Regression  analyses  of  the  m^  anomalies  vs  station  crustal 
amplification  factors  (After  Der,  McElfresh  and  Mrazek  1979) 
This  demonstrates  that  crustal  amplification  cannot  explain 

the  anomalies. 

Ms-mb  1'elat  ionships  for  NTS  and  Kazakh  at  AWRE  stations. 
(After  Marshall  and  Basham  1972). 

11  Results  of  some  theoretical  simulations  for  deriving  a  t*  vs 
m^  dependence.  The  m^  ■  1.35  t*  line  is  superposed. 

12  Map  of  stations  used  in  the  "NTS  experiment".  (NTSE) 

13  Averaged  relative  trace  amplitudes  derived  form  the  NTSE. 

14  Relative  t*  results  from  the  NTSE. 
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Station  t*  terms  for  selected  North  American  stations  derived 
trow  Soviet  nuclear  explosion  daLa  (After  Shore  1983). 


16  Patterns  of  S  waveforms  across  the  United  States.  The 

ch.-np.os  in  periods  and  amplitudes  indicate  low  Q  under  the 
southwestern  U . S  . 

17a-b  Typical  spectra  of  nuclear  explosions  from  NTS  and  Kazakh  as 

observed  at  NORSAR.  a)  Kasseri,  b)  10  Feb  1972  explosion  at 

Kazakh  . 

18  Results  of  regression  analyses  of  worldwide  t*  data  from 

nuclear  explosions.  Station  terms  are  on  the  left,  test  site 

terms  on  the  right. 
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19  Inconsistency  in  the  results  of  discrimination  analyses  for 
two  suspected  explosions  near  the  Baykal  rift  zone.  (after 
Savino  et  al  .  ,  1975)  . 

20  Body  wave  mangitude  anomalies  vs.  P  velocity  (After  Marshall 

et  al.,  1979).  n 

21  S-SS  travel  time  differences  for  the  U.  S.  (top)  and 
worldwide  (bottom)  after  Burdick  (1982). 
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Obscrved-to-souree  spectral  amplitude  ratio  shapes 
implied  by  the  proposed  versus  frequency  curves  and  the 
absorption  band  models  in  Fig.  Sv. 
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Station  effects  obtained  from  analysis  of  WW8SN  observations  together  with  effects  estimated 
by  Booth  el  al.  (1974)  for  the  LRSM  stations.  Contours  divide  areas  of  generally  positive  effects  from 

generally  negative  effects. 
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Magnitude  residuals  for  Booth  et  al  (1374)  plotted 
against  the  logarithms  oi  crustal  amplification  factor 
A  .  The  data  points  tent  to  cluster  around  two  regres 
s£on  lines,  one  for  the  :’.US,  the  other  for  the  WUS. 


Decrease  of  hl  computed  from  synthetic  pulses  for 
10  kt  and  a  100  kt  nuclear  explosion  due  to  incre; 
t*.  The  empirical  relationship  Am^  ~  1.35  At*  is 
sketched  in  for  comparison.  The  overall  slopes  o 
three  curves  are  similar. 
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Tracings  of  short-period  SH  phases  at  LRSM  stations 
across  the  United  States.  Depending  on  the  fre¬ 
quency  content  of  the  signals,  the  time  domain  mani¬ 
festations  of  anelastic  attenuation  vary,  but  the 
overwhelming  majority  of  the  signals  show  a 
diminution  of  amplitudes  and  the  decrease  of  high 
frequency  content  in  most  of  the  WUS  with  especially 
severe  effects  in  the  southwestern  United  States,  No 
corrections  for  radiation  patterns  were  made  in  these 
figures.  Instrument  gains  are  shown  on  each  trace. 
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Observed  relationship  between  Pn  and  log.^(A/T)  residuals, 
comparison  with  calculated  corrections  and  Qa  for  T  ■ 
(from  Marshall  et  al,  1979). 


l 

i 


.75s 


(a) 


•  Midpoint  in 
Tectonic  Region 

« 

*  Midpoint  in 
Shield  Region 

- - TNA 


* 

- —  — — x—  — -2—  .j?*—  £x-xx*>&^  — —  x— x  SNA 


5  50  55  6 


5  50  55  60 

Distance  (degrees) 

Plot  of  6t  (SS-S  time  residuals  in  seconds)  versus 
distance  for  (a)  North* American  data  and  (b)  foreign  test  site  data. 
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